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In this work we propose and analyse a scheme where the full spatio-temporal correlation of twin 
photons/beams generated by parametric down- conversion is detected by using its inverse process, 
i.e. sum frequency generation. Our main result is that, by imposing independently a temporal delay 
At and a transverse spatial shift Ax between two twin components of PDC light, the up-converted 
light intensity provides information on the correlation of the PDC light in the full spatio-temporal 
domain, and should enable the reconstruction of the peculiar X-shaped structure of the correlation 
predicted in [U|3] . Through both a semi-analytical and a numerical modeling of the proposed optical 
system, we analyse the feasibility of the experiment and identify the best conditions to implement 
it. In particular, the tolerance of the phase-sensitive measurement against the presence of dispersive 
elements, imperfect imaging conditions and possible misalignments of the two crystals is evaluated. 
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Introduction 

Recent theoretical investigations [T]-[3] outlined a pecu- 
liar spatio-temporal geometry of the biphoton correlation 
characterizing the entanglement of twin beams generated 
by parametric down- conversion (PDC). In collinear phase 
matching conditions, the biphotonic correlation displays 
a X-shaped geometry as a function of the relative spa- 
tial and temporal coordinates. This structure is non- 
factorable in space and time, offering thus the relevant 
possibility of manipulating the temporal bandwidth of 
the entanglement of twin photons by acting on their spa- 
tial degrees of freedom. The name "X-entanglement" was 
used pQ to describe this geometry. A key feature that 
emerged was the extreme spatial and temporal localiza- 
tion of the biphotonic correlation, on the micrometer and 
femtosecond range, respectively, which is present only 
when twin photons are detected in the near field of the 
PDC source. This feature allows in principle the gener- 
ation of ultra-broadband temporally entangled photons, 
via a proper control of their spatial degrees of freedom. 
We also showed that these features of X-entanglement 
persist in the high gain regime of PDC, where stimulated 
down-conversion becomes the main source of twin photon 
pairs. A detailed study can be found in [TJ [2] for type I 
phase-matching, and in [3] for type II phase-matching. 

The goal of the present work is a careful theoretical in- 
vestigation of a scheme based on the use of sum frequency 
generation (SFG) as a tool to explore the predicted X- 
shaped geometry of PDC entanglement. An experiment 
based on this scheme is currently under development at 
the Insubria University in Como [4] [5] . 

At low gains, a prominent way to probe the twin pho- 
ton correlation is the Hong-Ou- Mandel (HOM) detection 
scheme [6]. However, the experiment in Como works in 
the high gain regime of PDC, where the visibility of a 
HOM dip would be exceedingly low. We consider there- 
fore an alternative detection scheme based on the SFG 



process taking place in a second nonlinear crystal, onto 
which the PDC source is imaged. This second crystal op- 
erates as an ultrafast optical correlator, the up-converted 
field containing the information about the correlation of 
the injected source. 

Recent experimental works [7HTT] used SFG to test the 
twin photon/twin beam correlation in the purely tem- 
poral domain, by imposing a controlled temporal delay 
between the twin photons. A careful theoretical mod- 
eling of such schemes, and of related schemes based on 
two-photon absorption [12], has been developed in [TT] . 
However, this analysis is restricted to models considering 
only the temporal degrees of freedom of light. A fully 
spatio-temporal model for PDC/SFG has been consid- 
ered in [10], however this model is in turn valid only in 
the coincidence count regime of PDC. 

The aim of the experiment being developed in Como 
is to use SFG for exploring the PDC correlation in the 
whole spatio-temporal domain, by manipulating indepen- 
dently the temporal and the spatial degrees of freedom 
of the twin beams. To this end, we develop here a fully 
spatio-temporal model in order to describe the non linear 
processes taking place in the two crystals and the propa- 
gation between the crystals. This model is valid both in 
the high gain of PDC, where the experiment in Como is 
being performed, and in the low-gain regime, where twin 
photons can be resolved by coincidence counting. 

We demonstrate that the proposed scheme allows 
in principle the reconstruction of the X-shaped spatio- 
temporal correlation of twin beams/ twin photons, and 
we identify the best conditions under which an experi- 
mental observation can be performed. In particular, we 
analyse important issues such as the visibility of the in- 
formation thereby obtained (which becomes crucial in 
the high gain regime of PDC), and the tolerance of the 
scheme with respect to common experimental imperfec- 
tions, such as errors in the imaging scheme that maps 
the PDC light onto the SFG crystal or misalignments of 
the two nonlinear crystals. 
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Besides being necessary to describe the detection of the 
spatio-temporal X-correlation, we remark that our fully 
spatio-temporal model for PDC and SFG is also crucial 
to interpret the results of similar schemes aimed at ex- 
ploring the purely temporal correlation of twin photons, 
at least when the temporal bandwidths in play are large. 
This statement, further demonstrated by the experiment 
reported in [5], originates from the intrinsic nonfactora- 
bility of spatial and temporal degrees of freedom of the 
PDC correlation, which implies that the temporal prop- 
erties of twin photons depend also on the way their spa- 
tial degrees of freedom are manipulated before being up- 
converted (by the use e.g. of pupils that restrict t he ac - 
ceptance angle). A clear example is given in Sec. VII d, 
where we show that the effect of spatial diffraction in- 
troduced by an imperfect imaging scheme is a temporal 
broadening of the correlation, in addition to the more 
obvious spatial broadening. 

The paper is organized as follows. Sec Jl] illustrates the 
proposed scheme , while in Sec. [Qlwe develop a modeling 



of the setup. In Sec. [HI] we derive a general solution for 
the coherence function of the generated SFG light, so that 



in Sec IV we are able to demonstrate how the full spatio- 
temporal X-correlation can be retrieved by monitoring 
the SFG light intensity as a function of the temporal 
delay and the spatial transverse shift between the twin 
beams. Section VI discusses the issue of the visibility of 



the information, while Sec |VII| analyses the tolerance of 
the scheme with respect to some common experimental 
imperfections. 
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FIG. 1: Proposed scheme for detecting the full spatio- 
temporal correlation of PDC via sum frequency generation. 
The PDC light generated by a first crystal is imaged onto 
a second identical crystal, where up-conversion takes place. 
A transverse spatial displacement Ax and a temporal delay 
At are imposed independently on beam + and beam — re- 
spectively, by means of rotations/translations of the mirrors 
M1/M2 placed in the 2-f plane of the imaging system. The 
output SFG intensity is monitored as a a function of Ax and 
At 



I. GENERAL DESCRIPTION OF THE SCHEME 

The main features of the scheme we propose and of its 
theoretical modeling are the following: 

i) Both the spatial and the temporal degrees of free- 
dom of the optical fields are taken into account. In order 
to observe X-entanglement, a large temporal bandwidth 
(hundreds of nanometers) of PDC radiation needs to be 
collected [1, 2 , so that special care will be taken to model 
broadband field propagation. The spatial bandwidth is 
also rather large so that we assume that optical elements 
have large acceptance angles. 

ii) Our description is valid in any gain regime of PDC, 
so that it models both the generation and detection of 
twin-photons and of twin-beams. 

iii) A simplified scheme is illustrated in Figfl] A first 
nonlinear crystal, pumped by a broad coherent pump 
beam, generates PDC radiation. After eliminating the 
pump beam, the exit face of the PDC crystal is im- 
aged onto the entrance face of a second nonlinear crys- 
tal, where the inverse SFG process takes place. This 
is described in Fig{l] by a 4-f lens imaging system: in a 
real implementation, however, dispersive optical elements 
should be avoided as much as possible, since dispersion 
would drastically deteriorate the phase-sensitive correla- 
tion. The use of parabolic mirrors in place of lenses is a 



valid alternative, as we recently demonstrated in [HE]. 

iv) In order to have an efficient up-conversion process, 
the SFG and PDC crystals must be of the same material 
and cut for the same phase-matching conditions, while 
their lengths can differ. 

v) We focus on type I PDC, and in the following we 
shall always refer to the case of BBO (Beta-Barium- 
Borate) crystal pumped at 527nm in a e-oo phase- 
matching configuration. The PDC field is then described 
by a single field operator (ordinarily polarized) , while the 
pump field is extraordinarily polarized. 

vi) In order to explore the full spatio-temporal struc- 
ture of PDC entanglement, a key requirement is the abil- 
ity to impose independently a temporal delay At and a 
transverse spatial displacement Ax between the two twin 
components of PDC radiation. The twin components of 
the signal beams can be manipulated separately in the 
2f plane of the imaging device (the far-field plane with 
respect to both the PDC crystal exit face and the SFG 
crystal entrance face): the reason is that for a broad 
pump beam twin photons are always emitted with op- 
posite transverse wave- vectors q and —q, so that when 
a photon is found in the upper half of the 2-f plane, its 
twin will always be in the lower half plane. We shall de- 
note by beam + and beam -, respectively, the upper and 
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lower portions of the 2f plane, which correspond to the 
two portions of the PDC radiation having positive q y > 
and negative q y < components of the transverse wave 
vector along the y- axis in Fig{l] 

vii) Since the PDC correlation is strongly localized 
both in space and time (in the micrometer and fem- 
tosecond range respectively [IH3] ) , the temporal and spa- 
tial relative displacements of the twin beams must be 
scanned with micrometric precision. This can be realized 
by means of the two plane mirrors Ml and M2 placed in 
the 2-f plane of the telescopic system. As it will be de- 
scribed in detail in the following , a rotation of mirror 
Ml by an angle A<p generates at the SFG crystal input 
face a transverse displacement of beam + of an amount 
Ax = 2/A0, while a translation of mirror M2 by a dis- 
tance Az2 generates a temporal delay At = Az2/c. 

The intensity of the up-converted field generated by 
the SFG process in the second crystal is then monitored 
as a function of the relative temporal delay and spatial 
shift of the twin beams. The main idea that we will 
demonstrate is that this quantity is able to give a precise 
information about the structure of the correlation of twin 
beams generated in the first crystal and can be used in 
order to reconstruct the shape of this correlation in space 
and in time. 

In the treatment that follows we develop a model 
that takes fully into account the propagation effects and 
the phase-matching mechanism that selects the spatio- 
temporal frequencies in both crystals. 



II. MODELING THE OPTICAL SYSTEM 

In order to describe the scheme of Figjl] we consider 
separately the propagation in the PDC crystal (Step I), 
the linear propagation between the two crystals (Step II) 
and the up-conversion process in the second crystal (Step 
III). As indicated in Fig{l] the field operators in the dif- 
ferent planes of interest will be labeled with 

- aj at the input plane of the PDC crystal; 

- bj at the output plane of the PDC crystal; 

- Cj at the input plane of the SFG crystal; 

- dj at the output plane of the SFG crystal. 

The j = index refers to beams with central frequency 
ujq (either the pump field in the first crystal or the up- 
converted field in the second crystal), the j = 1 index 
to the beams of central frequency uj\ = uoq/2 (either the 
down-converted field in the first crystal or the fundamen- 
tal field in the SFG crystal). 

The description of field evolution along z, the mean 
propagation axis of the system, will be performed ei- 
ther in the direct spatio temporal space (5, £), where 
x = (x, y) is the 2D transverse coordinate and t is time, 
or in the Fourier spatio-temporal domain (g, Q) where q 
is the transverse component of the wave-vector and ft 
denotes the frequency offset from the central frequency. 
For convenience, we shall use a compact notation for 
these space-time coordinates by making the substitutions 



(q, Q) — > w 

q • x — iQt — »• w • £ 



(la) 
(lb) 
(lc) 



Our analysis will be carried out at two levels 
1) In the limit where the pump beam driving the PDC 
process is broad and long enough, we shall adopt the 
plane- wave pump approximation (PWPA), which allows 
us to derive analytical or semi-analytical results. This 
model will be presented in the next sections [Tl]|VT 



2) In order to obtain results for a finite pump, we also de- 
veloped a full 3D+1 numerical model, based on stochastic 
simulation of field evolutions. This model will be intro- 
duced in Sec. [VTl 



Step I: propagation in the PDC crystal 

In this section we describe the model we use to describe 
PDC (also derived in [TJ [2j [13]), and we recall the main 
features of the space-time correlation of PDC light that 
we called X-entanglement [IH3]. 

The pump and down- converted fields are described by 
two field operators bo and &i, centered around the fre- 
quencies ujq and uj\ = co>o/2, respectively. Normalization 
is such that (bj(£)bj(£)) gives the photon number per 
unit area and unit time. The generation of the PDC 
field along the nonlinear crystal takes its simplest form 
in the Fourier domain 



bj(w, z) 



(27r)i 



■bj(€,z)e 



j = o, l . 



(2) 



where we recall that w = q, f2 is the set of 3D Fourier 
coordinates, £ = (x, t), while z is the longitudinal coordi- 
nate along the mean propagation direction in the crystal. 
Next, we introduce the slowly varying amplitudes 



b {w,z) = e- iko ^ )z bo{w,z). 



(3a) 
(3b) 



where kj Z (w) = y/kj(w) 2 — q 2 is the z-component of the 
wave- vector for the j-field. These amplitudes vary slowly 
along the z-coordinate, because their evolution is only 
due to the nonlinear interaction, since we have subtracted 
the effect of the fast linear propagation contained in the 
phase factor exp[ikj Z (w)z]. We can assume the pump 
beam is undepleted by the PDC, so that its evolution 
is only linear bo(w,z) = bo(w,0) . In the same approx- 
imation, the pump operator can be substituted by its 
c- number amplitude bo(w,0) — >• a p (w). The propaga- 
tion equation for the signal field contains only first-order 
z-derivatives and takes the form [lj [13] 



9_ 

Yz 



h(w, z) 



-I 



xe 



dw 1 



(2tt) 3 / 2 

iA(w,w ') 



a p (w + w r )b\(w f < z) 



(4) 
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Here the phase matching function 

A(w,w') = k lz (w) + k lz (w') - k 0z (w f + w') (5) 

describes the phase-mismatch between the two generated 
signal modes (<f, Q), (<f and the pump mode (<f + 
<f', Q + Q'). Efficient down-conversion takes place only in 
those modes for which the phase-mismatch is small. The 
coupling constant a is defined by 



d. 



hir 3 u 



eff\ 



4eonon^c 3 



(6) 



where d e ff is the effective second order susceptibility of 
the nonlinear crystal, no and n\ are the refraction indexes 
at the central frequencies ujo and uo\. 

Eq.Q can be analytically solved in the plane- wave 
pump approximation, ap^w + w') — ^ (27r) 3 / 2 a p 5(w + w '), 
where a p denotes the pump field peak value in direct 
space. As analyzed in detail in [2], such an approxima- 
tion holds as long as the pump beam waist and dura- 
tion are larger than the spatial transverse displacement 
and temporal delay experienced by the pump and signal 
beams along the crystal because of walk-off and group 
velocity dispersion (in the example of a 4mm long BBO, 
a pump pulse with a waist larger than ~ 250 /im and 
a duration above ~ 200 fs satisfies this condition). The 
solution is expressed by a unitary transformation link- 
ing the field operators at the output face of the crystal, 
bi(w) = b\{w, z = / c ), to those at at the input face, 
a\[w) = bi(w, z = 0): 



bi(w) = U(w)ai(w) + V(w)a[(— w) , 



(7) 



The explicit expression of the functions U and V can be 
found e.g. in [2 . Here we notice that they depend on w 
only through the plane-wave phase mismatch 

A pdc (w) = A(ixT, -w) = k lz (w) + k lz (-w) - k (8) 

where ko = uquoq/c is the wave number of the pump. 
All the properties of the PDC light are described by the 
following second-order field correlation functions 



(&{(w)6i(w')} =S(w 
(b 1 (w)b 1 (w , )) =S(w 



w')\V(w)\ 2 (9a) 
w')U(w)V(-w) , (9b) 



In particular, from Eq. (9b), we see that the function 

Fp DC (w) = U(w)V(-w) (10) 

represents the probability amplitude that a pump pho- 
ton at (q = 0,0 = 0) is down-converted into a pair of 
phase-conjugated photons (<f, Q) and (—q,—Q). Its ex- 
plicit expression is 



^pdc(^) ge 



ikol, 



sinhr(«J)^ c 



T(w)l, 



|coshr(ty)l c 



2T(w)l c 



sinhr(it>)Z c 



T(w)l c 



[AP dc (^)y 2 



(lla) 
(lib) 



where g is the dimensionless gain parameter proportional 
to the pump peak amplitude 



9 



rl r a rj 



(12) 



The other relevant function is the PDC spatio-temporal 
spectrum 



7?\|2 



\V(w)\ 



2 sinh 2 [T(w)l c 
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(13) 



which gives the photon number distribution in the spatio- 
temporal Fourier domain. 

We remark that equations p|-([l3|) are valid in any gain 
regime of PDC. In the low gain regime ^ C 1 they de- 
scribe the down- conversion of pump photons into pairs 
of signal-idler photons that can be resolved individu- 
ally. In the high gain regime g ~ 1 stimulated down- 
conversion becomes important, and these equations de- 
scribe the generation of macroscopic twin beams of light 
made of "bunched" pairs of twin photons. 
In particular, in the low-gain regime Eq. (11) takes the 



well known Sine dependence on the phase mismatch: 



limFp D c(^) = ge ik ° lc Smc 



A pdc (w)l c 



(14) 



In Ref. pQ the name X-entanglement was used to de- 
scribe the shape of the spatio-temporal correlation of the 
biphoton amplitude at the crystal output face. The quan- 
tity of interest is therefore: 



^dc(0 = <MO&i(£ + ao> 



(15) 



which in the stationary plane- wave pump regime depends 
only on the relative spatial and temporal coordinates 
A£ ee (Ax, At). It can be expressed [2 as the Fourier 
transform of the spectral probability amplitude Fpdc(^) 
of generating photons in phase conjugate modes w and 
—w , that is 



^PDC(0 



/ 



dw 



(16) 



Fig. [2] shows the behaviour of |Fpdc| fc> r collinear phase 
matching, calculated by using the Sellmaier relations for 
the refractive indexes [16 j. As can be inferred from Eqs. 
( pT| and ( [l4| ), in any gain regime Fpdc(^) is strongly 
peaked around the phase matching curves defined by 
A pdc (w) = 0. The hyperbolic geometry displayed in the 
neighborhood of degeneracy can be understood by mak- 
ing a quadratic expansion of k z around (q = 0, Q = 0) 



k lz (q, + k[n + -k"Vt 2 



q 

2fci 



(17) 



where k\ = k\ (w = 0), k[ = dki/dQ\^ = o and k'{ = 
d 2 kj / dVL 2 \yt=o • 1 The phase- matching function (M) takes 



1 A detailed discussion on the range of validity of the approximations 
( |17| > can be found in 3 in the context of type II phase-matching. 
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FIG. 2: Plot of IFpdcI in the (q x , Q)-plane, for a 4mm long 
type I BBO pumped at 527nm, cut for collinear phase match- 
ing. The unshaded region corresponds to the frequency filter 
used in the simulation with 2Q max = 0.965 x 10 15 Hz FWHM. 
The band widths ~ 4qo and ~ of the phase-matching of 
I-FpdcI along the spatial and the temporal frequency axis are 
indicated. The parametric gain is g — 8, Aq DC = 0. 
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then the quadratic form 



ql U? D 



(18) 



where Aq DC = 2k\ — ko is the collinear phase-mismatch 
parameter, and 



D 



Qd 



\AiAc 



(19) 



qjj and Qjj determine the characteristic scale of varia- 
tions of Fpdc along the ^-axis, at fixed and along 
the g-axis at fixed respectively. They scale with 
the inverse square root of the crystal length and are 
generally much smaller than the range of frequencies 
of the whole PDC emission spectrum (the latter can 
in principle extend up to the pump optical frequency). 
When tuning the crystal for collinear phase matching, 
i.e for Aq DC = 0, phase-matching occurs along the lines 
q/qd — dzfi/n^). Under this condition, we no tice t hat the 
first zeros of the function T(q,Q)l c (see Eq.(llb)) along 
the Q— and q— axis, evaluated using approximation (18) 



& D = 



V^ 2 V) 1/4 ^D, to = V / 2(^ 2 V) 1/4 to, (20) 

provide good estimates of the widths of |Fpdc| along 
those axis, which takes into account the gain broaden- 
ing effect in the spectral domain. For the example shown 
in Figj|J with g = 8, we have tip ~ 4Qd, qp « 4^, as 
indicated by the arrows. 

The hyperbolic geometry of phase matching in the 
spectral domain turns into a characteristic X-shaped 
geometry for the spatio-temporal biphoton correlation 
^pdc(Ax, At) ( [16] ) . Figure^ shows the profile of ^pdc 
as a function of the relative spatial and temporal co- 
ordinates, (only one transverse dimension is shown). 
The tails of the structure are oriented along the lines 



FIG. 3: Example of X-entanglement for the same crystal as in 
Fig. [2] (a) Plot of the modulus of the biphoton amplitude as 
a function of the relative spatio-temporal coordinates Ax, At. 
Temporal (b) and spatial (c) profile of the central peak of the 
correlation. Same parameters as in Fig [2] 



QpAt = ±qpAx : a feature expressing a linear relation 
between the temporal delay and the spatial transverse 
separation acquired by the twin photons when arriving 
at the crystal output face [UG], i.e. making use of the 
equations (19) 



At = ± y / k / {k 1 Ax 



(21) 



As described in detail in [2], when tilting the crystal away 
from collinear phase matching, this linear relation be- 
comes less and less stringent: as a consequence the tails 
of the structure become progressively less visible. The 
Aq DC = phase-matching condition represents therefore 
the optimal configuration to observe the X-entanglement. 

A relevant feature of the X-entanglement is the strong 
localization of the central peak of the correlation func- 
tion, whose cross-sections along the temporal and the 
spatial axes are shown in the lower part of Fig (3] . 
The width of the peak is in principle determined by 
the inverse of the full PDC emission bandwidth, or in 
practice by the bandwidth intercepted in the measure- 
ment. In the example of the figure, we simulated the 
presence of a filter in the temporal frequency of width 
2A^ max = 0.965 x 10 15 Hz (indicated by the unshaded 
region in FigJ5|, which corresponds to a wavelength inter- 
val ranging from 830nm up to the conjugate wavelength 
1444nm. The temporal profile of the correlation in Fig(3] 
has a full width half-maximum Atf w h m ~ 5fs, close to 
that of the Fourier transform of a box function of width 
2^max- The spatial width of the correlation peak, in turn, 
is determined by the range of spatial frequencies involved 
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in the PDC emission, and typically is on the order of a 
few microns. 



Step II: propagation between the two crystals 

We assume ideally that the optical setup illustrated in 
Fig |l] behaves as a perfect imaging system, free from dis- 
persion and losses. In the absence of any temporal delay 
and spatial shift (At = 0, Ax = 0), the PDC field at the 
PDC crystal exit face is then mapped into the SFG crys- 
tal entrance face c\ (x, t) = b\ (x, t) 2 . In order to model 
the effect of a small rotation of mirror M\ and a small 
translation of mirror M2 we assume that the propagation 
angles in play are such that the paraxial approximation 
holds. 

Let us first consider the temporal delay applied to the 
q y < Fourier modes of the PDC light (beam -). By 
translating the mirror M2 by Az2 the separation between 
the two lenses in Fig. [I] becomes 2f + Az2. The fields 
at frequency uji + Q at the left and right focal planes 
of the imaging device are related through the algebraic 
transformation (the x — »• — x reflection is omitted) 



ci(x, Q) 



f 2 



Az 2 



h&n) (22) 



The effect of diffraction due to the additional propagation 
Az2 is described, within the paraxial approximation, by 
the quadratic phase factor proportional to \x\ 2 /f 2 . Its 
effect can be neglected as long as 



Az < 



A f 



7T \X 



12" 



(23) 



Referring to the conditions of the experiment [H [5] , for 
/ ~ 20cm, A ~ If mi and \x\ ~ w p < 1mm the condition 
reads Az <C 1.2cm. In order to explore the X-shaped 
PDC correlation, delays of a few hundreds femtoseconds 
at most are sufficient (see e.g. the plot in Fig. [3^i) , so that 
in practice condition (23) is always fulfilled. Each Fourier 
mode of beam - undergoes therefore, within a very good 
approximation, the following diffractionless transforma- 
tion 



ci(g,fi) 



l Az 2 



bi(q, Q) for q y < 



(24) 



where an inessential constant phase factor has been omit- 
ted. 

We now consider the manipulation of beam + through 
the rotating mirror Ml. As shown in AppjAj a rotation 
of Ml by a small angle Acj) around a given axis gener- 
ates a transverse shift As = 2/A0 of the beam at the 
entrance face of the SFG crystal (the imaging plane), in 



2 For simplicity we neglect here and in the following the spatial re- 
flection with respect to the z-axis x — > —x, as well as constant 
phase factors 



the direction orthogonal to the rotation axis. According 
to Eq. (A6), the complete transformation undergone by 
beam + between the input and output planes of the 4-/ 
telescopic system can be written in the Fourier domain 
in the form (the q —> —q reflection and the minus sign 
are here omitted for simplicity) 



d(q, n) = e^ A %(q, SI) for q y > 



(25) 



where As denotes the transverse shift generated by a 
rotation of mirror Ml around a generic axis. Putting to- 
gether relation (24) and (25), the overall transformation 



describing propagation from the PDC crystal output face 
to the SFG crystal input face can therefore be synthesized 
into the following unitary input-output relation 

Cl (q, SI) = [Q(q y )H + (q, SI) + <~>(-q y )H_(q, SI)} h(q, SI) 

(26) 

where &(q y ) denotes the step function, equal to 1 for 
q y > 0, to zero for q y < 0, and 

H+(q,n) = e^ As ; H_(q,n) = e^ Az * (27) 

are the transfer functions associated to beam + and beam 
— , respectively. The identity Q(q y ) + 0(— q y ) = 1 guar- 
antees that the commutation rules are preserved by the 
transformation. 

It is important to notice that the correlation measure- 
ment depends on the particular choice of the rotation axis 
of Ml, as it will be shown in SeclIVI We shall consider 



explicitly the two configurations illustrated in Fig(4j in 
case a) Ml is rotated around the 7r-axis orthogonal to the 
gap between the two mirrors, while in case b) the rotation 
axis 7r of Ml is coincident with the x-axis (orthogonal to 
the figure plane). We have thus 



As = (Ax,0), Ax 
As = (0,Aj/), Ay - 



■ 2fA(j) x in case (a) (28a) 
2fA(/) y in case (b) (28b) 



where A(j) x and A(j) y denote a small rotation angle ap- 
plied to the mirror around the 7r-axis in configuration (a) , 
around x-axis in configuration (b). 



Step III: propagation in the SFG crystal 

We now model the generation of the up-converted field 
inside the SFG crystal. This process is just the reverse of 
the down- conversion described in Sec. [TTJ Following the 
same procedure outlined there, we work in the Fourier 
domain and we introduce the slowly varying amplitudes 
of the fundamental d\ (carrier frequency cji), and second 
harmonic do (carrier frequency c^o) fields along the SFG 
crystal 

do(w,z) = e- ikoz ^ )z d (w,z) (29a) 
d ± (w,z) = e- ikl 'W z d!(w,z), (29b) 
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q(5,n) 



Equation (30a) describes all the up- conversion processes 
where a pair of fundamental photons in modes w f , w — w' 
are up-converted into a SFG photon in mode w. Ac- 
cordingly, in this equation A(w',w — w f ) — k\ z (w r ) + 
k\ z {w — w') — ko z (w) represents the phase mismatch of 



such a process in the second SFG crystal. Equation ([30b]) 
is obviously just the same as Eq. Q describing the down- 
conversion process. 



Ml 



M2 



Equations (30) need to be considered together with 
the initial conditions at the crystal input face. For the 
fundamental field we have 



(b) 



SFG cr\stiil 



FIG. 4: In configuration (a) the rotation by a small angle 
A<j) x of mirror Ml around the 7r-axis produces the displace- 
ment As = (Ax = 2/A0 x ,O) along the x-direction at the 
entrance face of the SFG crystal (orthogonal to the figure 
plane). In case (b) the rotation is performed around the x- 
axis and produces the displacement As = (0, Ay = 2fA<j) y ) 
along the ^-direction. 



di(w 1 z = 0) = Ci(w) (31) 
where c\ is the field down- converted in the first PDC 
crystal, after Step I and II. For the second harmonic field, 
we assume that the pump field is completely eliminated 
after the PDC crystal, so that it is in the vacuum state 
at the SFG crystal input. 



which vary along z only because of the nonlinear interac- 
tion, since we have subtracted the effect of linear propa- 
gation. Their evolution is described by the following pair 
of coupled equations: 



ddo(w, z) 
dz 



ddi(w, z) 
dz 



= —a 



w — w ' , z) 



xe 



iA(w ,w — w )z 



(30a) 



r dw' 



do(w + w\ z)d\(w\ z) 



)3/2 



xe 



iA(w,w ')z 



(30b) 



We now assume that the SFG crystal is short enough 
so that only a small fraction of the PDC light is up- 
converted. In these conditions, the fundamental field 
remains basically unchanged during propagation, and 
Eq. ( 30 ) can be solved with a perturbative approach sim- 



ilar to that used in the low gain regime of PDC in Ref.[2] 
(see App.A therein). In this way, we obtain an explicit 
expression that links the operators at the SFG crystal 
output plane to those at the SFG input plane: 



di(tZT) = d 1 (w,l , c ) = e iklz ^c 1 (w), 



do(w) = d (w,l f c ) = e ikoz ^ 



co(w) 



dw 1 



(2^) 3 / 2 



ci(w — w ' ) C\ (w r ) Fsfg ( w — w' ,w') 



(32a) 
(32b) 



where The function Fsfg (w, w f ) can be interpreted as the prob- 

FsfgK^O = <e i A( ^ /) ^ Sinc A( ^ ( f ' )l ' c . (33) 
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ability amplitude density 3 that a pair of photons in the 
fundamental modes w = (g, O) and w' = (<f , fV) are 
up-converted into the second-harmonic mode + w f = 
(q + (f , Q + Q'): this up-conversion probability is non neg- 
ligible only for those pair of modes for which the phase- 
mismatch A(w,w f )l f c « 0. As it can be expected, the 
probability amplitude for such a process is formally iden- 
tical to that of the reverse process of down- conversion. 
Eq.(33) has indeed the same form as Eq. (14), the main 



Inserting the input output relation (32b) for the SFG 



difference being that Eq. (14) describes PDC only in the 
plane wave pump limit, where the only allowed down- 
conversion processes are those leading to twin photons in 
modes w and w' = —w. 



III. GENERAL SOLUTION: THE SFG 
COHERENCE FUNCTION 

We now put together the chain of field transformations 
presented in the previous section. Our main goal will be 
to evaluate the intensity of the SFG field, but we start 
from a more general result, i.e. the coherence function 
of the SFG field in the spectral domain, evaluated at the 
output face of the SFG crystal 

Csfg(^o,^o) = (4(™o)do(w' )) • (34) 

3 Its square modulus \Fqfq(w,w f )\ 2 gives the probability per unit 
of spectral bandwidth. 



crystal inside this expression, we obtain an equation that 
links the SFG coherence function to the correlation func- 
tions of the fundamental field c\ at crystal input face. 
These correlations can be calculated by using relations 
^ together with the identities \H+(w)\ 2 = \H_(w)\ 2 = 1 
and @(q y ) + 0(— q y ) = 1. We obtain the following ex- 
pression 



(c{(w)c 1 (w f )) = 5(w-w')\V(w)\ 2 (35a) 

(c 1 (w)c 1 (w , )) = 6(w + w')U(w)V(-w) 

x [Q(q y )H + (w)H_(-w) 

+ e(-q y )H+(-w)H-(w)] (35b) 



Using these relations we obtain after some manipulations 
the following result 



Csfg(^o,^o) = C^(w ,w' ) +4fg K,^o) (36) 



with 



C^(w ,w' ) = 2S(w - w' ) I j^lViwo - w)\ 2 \V(w)\ 2 \F SFG (w - w,w)\ 2 



(2tt) 3 1 
dw 



(2tt) 3 / 2 



G(q v )H + (w)H_(-w)Fp D c(w)F SFG (w, -w) 



(37a) 
(37b) 



where Fpdc(^) and Fsfg(w, w ') are the spectral prob- 
ability amplitude defined in E qs.([TT|) and (33). 

The incoherent term ( |37a| ) originates from the up- 
conversion of PDC photons which do not belong to phase- 
conjugated mode pairs. By contrast, the coherent term 
(37b) originates from the up-conversion of photon pairs 
coming from phase-conjugate PDC modes, this latter 
process leading to the partial reconstruction of the origi- 
nal coherent pump beam as described in the experiment 
in Ref.[4 . In this regard, it is convenient to introduce a 
special symbol for the up- conversion probability ampli- 
tude appearing in Eq.( |37b| ), namely 



^sfg(^; -w) 



alt 



(38a) 

Sinc^« , (38b) 



where 



A sfg (w) A(w, -w) = ki(w) + ki(-w) - k { 



sfg 



(39) 



is the phase mismatch function for these coherent up- 
conversion processes. The function F^™^ plays a fun- 
damental role in determining the properties of the co- 
herent component of the SFG field. It describes the ef- 
ficiency with which photons belonging to a particular 
pair of phase-conjugated modes (w) and (—w) under- 
goes the inverse process of the original PDC event, i.e. 
the back-conversion into the monochromatic plane-wave 
mode (w = 0) corresponding to the original pump driv- 
ing the PDC. The corresponding phase- mismatch A sfg is 
identical to the PDC phase-matching function ^ except 
for the presence of the wavenumber &Q fg , which differs 
from the pump wave number ko only in the case when the 
two crystals are not perfectly aligned (for the type I e-oo 
phase-matching we are considering only the pump beam 
k- vector depends on direction). When the two crystals 

are aligned, i^FG^ coincides with the probability ampli- 
tude for the reverse PDC process in the low-gain regime, 
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as given by Eq.(|14|) (apart from a con stant factor). 

As can be inferred by comparing Eq.( |37b| ) and Eq.(16), 
only the coherent contribution contains the information 
about the biphoton correlation ^pdc(0 we are looking 
for. The incoherent contribution rather acts as a back- 
ground which tends to deteriorates the visibility of the 
correlation measurement HTl. 



Notice also that the coherent component (37b) factor- 



izes into the product of the SFG field mean values: 



C^\w,w') = (4(w))(d (w')) 



(40) 



with 



(d (w)} = S(w) 



dw' 



-H + (w')H_(-w') (41) 



(2tt) 3 / 2 



From Eqs.([36|)-(37) we can evaluate the SFG intensity 



(number of photons per unit area and unit time) at the 
exit face of the 2nd crystal, obtaining 



wo = (4(i)dod)) 



r(coh) 
SFG 



r (inc) 
Z SFG 



(42) 



with 



r (inc) 
SFG 



r (coh) _ 
SFG — 



dw 
(2^)3 

dw' 



dw' 
(VP 



A|2 



\V(w)\'\V(w')\'\F SFG (w,w')\ 



Q(q y )H + (w')H_(-w')F PDC (w')F^\w') 



(43a) 
(43b) 



The intensity distribution of the SFG light in the near- 
field is uniform in space and time, an artifact due the 
spatio-temporal invariance of our model deriving from 
the monochromatic and plan e-wave p ump approxima- 
tion. By comparing Eq.(37b) and Eq.(43b), we obtain 



the following relevant identity 



C^(w,w') = (2nfS(w)S(w')I^ 



(44) 



which shows that the coherent component of the SFG 
field contains the same kind of information on the bipho- 
ton amplitude both in the space-time domain and in the 
spectral domain. The most important difference lies in 
that the coherent intensity in the spectral domain is con- 
centrated in a single peak at the origin (q = 0, O = 0), 
i.e. in the mode corresponding to the original pump 
field. On the other hand, the spectrum of the incoherent 
background is delta-correlated in space and time, and 
spreads over a very broad range of spatial and tempo- 
ral frequencies, its particular shape being related to the 
phase-matching conditions inside the SFG crystal as it 
will be elucidated in [14 . This circumstance suggests 
that either far-field or spectral measurement may be con- 
veniently used to enhance the visibility of the coherent 
contribution with respect to the incoherent background. 



IV. RETRIEVING THE PDC CORRELATION 

Let us now investigate how the information on the cor- 
relation ippr>c(%,t) °f twin beams/ twin photons can be 
effectively extracted fro m th e coherent component of the 
SFG light, given by Eq.(43b). Making explicit the depen- 



dence on the temporal delay and spatial shift applied to 
beam + and beam - contained in the transfer functions 



(27), it can be written in the form 



r (coh) 
Z SFG 



[As, At] = |VWs (As, At) f 



where 



^mea S (A|) =2 



/ 



dw 



e(q v )e M ^F PDC (w)F, 



h) 

SFG 



(45) 
(46) 

and Af (As, At). 

In this section we consider the configuration shown 
in Figj4^i, where mirror Ml is rotated orthogonally to 
the gap between the two movable mirrors so that As = 
(Ax = 2fAcj) x ,0). In this case, the r.h.s. of Eq.(45) coin- 



cides with the square modulus of the following quantity 



VWs(Ax,0,At) 

dn 



(47) 



J (2^)2 J 2^ 6 F PBC (qAl)F SFG (q,tt). 



The step function &(q y ) could be eliminated from Eq.(46) 
by exploiting the parity of integrand with respect to q x , 
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(a) 



, / s f G [o,o,o] 

oUtf ^-100 

At(fs)° -50 100 Ax( ^ m) 




At(fs) 



represents the probability amplitude that a pair of pho- 
tons delayed by a time At and displaced by As are co- 
herently back-converted into the plane-wave pump mode 
(g = 0, Cl = 0). This function can be also interpreted as 
the optical response function of the SFG crystal in the 
measurement of |^pdc| 2 via SFG. 



A. Thin SFG crystal limit 



FIG. 5: (Color online) (a) Surface and (b) density plots 
Jsfg[Ax,0, At] = |^meas(Ax,0, At)| 2 normalized to its peak 
value, corresponding to measurement in the configuration of 



the latter being a consequence of the radial symmetry 
of .FpDc(<?> 0) an d Fgp^(g, in the (q x , g^-plane. The 
function ^ meas (Ax, 0, At), like ^pdc, is therefore an even 
function of both the spatial and the temporal coordi- 
nates, a feature not holding for the alternative configu- 
ration with mirror Ml rotating around the x-axis (see 
FigHp) which will be discussed in the next section. Fig- 
ure]^] plots the coherent contribution of the SFG inten- 
sity as a function of the applied spatial and the temporal 
shifts. In this example, the two crystals have the same 
lengths l c = V c = 4 mm and the PDC gain is g = 8. 
The structure preserves the nonfactorable X-shaped ge- 
ometry of the original PDC correlation function |^pdc| 2 
shown in Fig(3j with symmetric tails developing along 
the bisector lines At = ±^k"kiAx (corresponding to 
the diagonals of the plot frame). 

The SFG coherent component (|47| ) differs from the 
PDC biphoton correlation I^pdc I ^defined in Eq.(16) 



only for the presence of the SFG spectral amplitude 
F^q\w). The latter describes how the phase match- 
ing mechanism selects the spatio-temporal modes in the 
coherent up-conversion process. We shall analyse in de- 
tail how it affects the PDC correlation measurement in 
the next subsection. In this regard, it is useful to recast 
Eq.(46) in the form of a convolution of the PDC bipho- 
ton amplitude ^pdc - the quantity under investigation - 
with the Fourier transform of Q(q y )F^ > Q > (q^Q)^ i.e. 

V'mea S (AO = f *fVsF G (A£ - f )V>PD C (f )• (48) 

where 



$3FG(A0 



/ 



(2tt) 3 



(49) 



This identity has been written taking into account that 
As = (Ax,0), so that the step function could be elimi- 
nated as for Eq.(47). Noticing that Eqs.(49)-(38) are for- 



mally identical to Eqs.( 16)-( 14), we see that ^sfg(A£) 
is proportional to the low gain biphoton amplitude in 
the space-time domain for the second SFG crystals. It 



Let us assume that the two crystals are tuned for 
the same phase matching conditions, so that A sfg (u>) = 
A pdc (uT), but the SFG crystal is much shorter than the 
PDC crystal V c « l c . In this limit, the PDC biphoton 
amplitude can be exactly reconstructed by monitoring 
the SFG coherent component. 

This can be shown by inspection of Eq. ( [46] ) . When the 
same phase matching conditions hold in the two crys- 
tals, the spectral probability amplitudes F$™q\w) and 
Fpdc(w) are peaked around the same geometrical curve 
A^ dc (w) = A s ^(w) = 0. However, for l' c « l c the 
spectral bandwidths in the SFG crystal are much wider 
than in the PDC crystal (the bandwidths of F$fg ex- 
ceeds those of Fpdc by a factor \fhJV c ^> 1, according 
to Eq.(19). As a result, F$™q\w) is almost constant in 



the region where Fpdc(^) is not negligible, close to its 
maximum value <jV c . Under these conditions, the filtering 
effect due to phase-matching in the SFG crystal becomes 
ineffective and we have 

/^[Aar.0, Az] « a% 2 |Vpdc(Ax, 0, At)| 2 (l' c « Z c }50) 

Thus the coherent component of the SFG output repro- 
duces the biphoton correlation as anticipated. 
The same conclusion can be derived also looking at Eq. 
(|48|): in the limits << Z c , the biphoton amplitude i^sfGi 



defined by Eq. ( 49 ) , has scales of variation in space and 
time much shorter than ^pdc? because it is the Fourier 
transform of a much wider spectral amplitude. It behaves 
therefore as a (5-function inside the convolution integral 
(48), so that the PDC biphoton amplitude is recovered. 



For completeness, we men tion that in the V c — » limit 
the incoherent contribution ( |43a[ ) takes the form 

ff 2 £|(&J(£,t)&iOM)>l 2 (I'c^o) (5i) 



r (inc) 
SFG 



where 



(blix^b^t)) = Jj^J (f^(^)l 2 (52) 

is the PDC photon flux evaluated at the output face of 
the first crystal [as can be easily inferred from relation 
(9a)]. However, we shall see in Sec VI that the temporal 
walk-off between the PDC field and the generated inco- 
herent SFG field plays a fundamental role in the forma- 
tion of the incoherent component even when considering 
short propagation distances. In practice, the validity of 
expressions (51) fails as soon as the finite length of the 



SFG crystal is taken into account. 
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(a) / c '=100|Lim 



temporal cross-section 
for Ax=25jum 



j[A^,Ar = 30fs]} 




i 

0.5 



At(fc) v 50 -100 Ax( ^ m) 

(b) / c =2000|Lim 




0.06 -i 
0.04 
0.02 
0.00 : 

0.06 
0.04 
0.02 



-75 -50 -25 25 50 75 
At (fs) 



-50 ^^^^^F" 

At(fs)° 50TlO0 Ax(iim) °- 00 




(c) / c =4000|Lim 



75 -50 -25 25 50 75 
At (fs) 




-50 ^^^^^ 

At (ft) 50 -100 Ax( ^ m) °- 00 . 



•75 -50 -25 25 50 75 
At (fs) 



FIG. 6: Left panel: numerical evaluation from Eq.(45)-(47) of 
/sfg[Ax, 0, At]/Jg& [0, 0, 0] for increasing values of the SFG 
crystal length l' c . Right panel: corresponding temporal cross 
section /sfg[Ax = 25/xra, 0, At]/I§§k[0, 0, 0]. The X-shaped 
structure is robust against propagation in the SFG crystal 
and the tails are enhanced with respect to the central peak 
as l' c increases. The PDC parameters are g = 8, l c = 4mm. 



B. Long SFG crystal 

We now investigate how propagation in the SFG crys- 
tal affects the retrieval of the PDC correlation via the 
measurement of [Ax, 0, Az], showing that both the 
X-shaped geometry and the strong localization of the 
biphoton correlation are preserved, at least when the two 
crystals are equally tuned. 

The left panels of Fig(6] plot the retrieved quantity 
I^fq [Ax, 0, At] for increasing values of the SFG crys- 
tal length V c . The most evident effect is that the tails of 
the correlation becomes progressively more visible with 
respect to the central peak, as the propagation distance 
inside the SFG crystal increases. This is clearly shown by 
the temporal cross-sections at Ax = 25/im of the same 
quantity (right panels of Fig{6]). Figure[7] shows that the 
height of the tails, at a fixed time delay At — 30 fs, 
increases almost linearly with V c up to the value of the 
PDC crystal length l c = 4 mm. We expect therefore that 
the choice of a few millimeter long SFG crystal, with re- 
spect to that of a very short crystal, presents a twofold 



0.03 
0.02 
0.01 
0.00 



OT,0] 



2000 4000 6000 8000 
V (urn) 



FIG. 7: Plot of the height of the tails of the retrieved corre- 
lation maxAcc {/sfg[Ax, 0, At = 30fs]}, at a fixed time delay 
At = 30fs, normalized to the peak value /sfg[0 5 0? 0] 5 as a 
function of l' c . 



advantage: i) the total number of up-converted photons 
is obviously larger and ii) the visibility of the tails, al- 
though small with respect to that of the central peak, will 
be enhanced, a feature which should facilitate the exper- 
imental observation of the X-shaped structure through 
the scanning of Ax and At. A further advantage lies 
in that the overall visibility of the coherent component 
with respect to the incoherent background improves sub- 
stantially as the SFG crystal length is raised above a few 
hundred of microns, as it will be shown in SeclVI 



An explanation of the behaviour depicted in Fig(6] 
can be obtained considering the convolution integral in 
Eq.(48 ) and referring to the density plots of Fig. [8j which 
illustrate the effect of the convolution (48) under three 
different conditions: 

(a) V c <C l c ' the SFG crystal response function ^sfg 
behaves as a Dirac (5-function both in the spatial and 
temporal domain (see plot a2), so that the 2nd crystal 
works as an ultrafast correlator with instantaneous and 
localized response. The coherent component of the SFG 
intensity (plot a3) provides therefore the PDC biphoton 
amplitude square modulus I^pdcI 2 , in agreement with 
the short crystal limit result (50). 

(b) l c ~ \' c : if the SFG crystal has a length comparable 
to that of the PDC crystal, it displays a nonlocal spatio- 
temporal response. However, provided the two crystals 
are tuned for the same phase matching conditions, the 
response of the SFG crystal has the same geometrical X- 
shape as the original biphoton correlation [see plots (bl) 
and (b2)]. This does not come unexpected: as a mat- 
ter of fact, the function ^sfg is identical to ^pdc in the 
low PDC gain limit (apart from a constant factor). In 
the high gain regime of PDC considered in this example, 
the PDC biphoton amplitude displays a faster decay of 
the tails than its low gain counterpart, 4 but nevertheless 
keeps the same geometrical shape. Therefore ^ meas is the 



4 The faster decay of the PDC correlation V'PDC along its tails for 
high gain is due to the broadening of the spectral amplitude with 
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FIG. 8: (Color online) The left and central panels shows the 
density plots of |^ PDC (Ax, 0, At)| and |^ SFG (Ax, 0, At)\, re- 
spectively, considering either short or long crystals (the values 
of l c and l' c are indicated). The square root of the measured 
quantity, y/ Is°fg ~ |V ; meas(Ax, 0, At) | , obtained by convolv- 
ing ^pdc and ipsFG, is shown in the right panel. See text for 
a discussion. 



convolution of two functions which have almost the same 
X-shaped structure in the space-time domain, and the 
main effect of this convolution is to enhance the weight 
of the tails with respect to the central peak, as can be 
seen from plot (b3). Most-importantly from our point of 
view, the non-factorable X-shaped geometry of the PDC 
biphoton amplitude is preserved, as well as its strong lo- 
calization in time (or in space) when particular temporal 
(or spatial) cross-sections are considered, 
(c) l c <C l' c . Because of the symmetrical role played by 
^pdc an d ^sfg in Eq.(48), a configuration with l c <C V c 
would provide a direct measure of |^sfg| 2 , which is pro- 



portional to the PWP biphoton amplitude in the SFG 
crystal at low parametric gains (accordingly, the density 
plots (c2) and (c3) are nearly identical, as for the plots 
(al) and (a3) in case a). However, in this case the SFG 
photon flux would be strongly reduced because of the 
lower gain of the PDC field associated with short crys- 
tals. 



V. ALTERNATIVE DETECTION SCHEME 

We now focus on the alternative configuration (b) for 
mirror Ml shown in Figj4]3, considering rotations per- 
formed around the x-axis. In contrast to configuration 
(a) treated in the previous section, the step function 



Q(q y ) cannot be eliminated from Eq.(46) and the SFG 



coherent component is now given by the square modulus 
of 



^eas(0,A,,At)=2/^/^ 

with Ay 2fA(j) y . 



e iq y Ay - iQAt G(q y ) 



(53) 
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FIG. 9: (Color online) (a) Surface and (b) density plots 
7§fg[0 5 Ay, At] = |V>meas(0, Ay, At)| 2 normalized to its peak 
value, corresponding a measurement in the configuration of 
Figgji. 

Because of the presence of the step function, the even 
symmetry with respect to the spatial coordinate is lost. 
As shown in Fig|9j the retrieved correlation function in 
the (At, Ay)-plane display a V-shaped geometry: the 
tails extends along the lines At = ±y/k*[kiAy in the 
Ay > half-plane, while they disappears for Ay < 0. 
Moreover, compared to the result obtained in configura- 
tion (a) shown in FigjEJ the tails are strongly enhanced 
with respect to the central peak. From simple geometri- 
cal optics considerations, it can be shown that negative 
values of the rotation angle A<j) y of mirror Ml around the 
x-axis prevents the overlap of beam + and beam - inside 
the SFG crystal, so that the up-conversion efficiency be- 
comes much lower under this condition (i.e for Ay < 0). 
On the contrary, a negative rotation A(j) y increases the 
overlap of beam + and beam -, enhancing thereby the 
up-conversion efficiency compared to configuration (a) in 
the Ay > region. 

Noticing that the complex amplitude ^ meas (Ax, 0, At) 
and ^ me as(0 7 Ay, At) given by Eq.(47) and Eq.([53| are 
linked through the relation 



5 (Ax,0,At) = -[^ 



(0, Ay = Ax, At) 
(0, Ay = -Ax, At)], (54) 



g described by Eq.|20} 



from which it can be inferred that the temporal profile 
across the central peak [obtained by setting Ax = in 
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FIG. 10: (Color online) Temporal cross-sections of the SFG 
coherent component obtained by varying At for (a) As = 0, 
(b) As = 50/xm and (c) As = — 50/xm. The solid line cor- 
responds to configuration (a), the dotted dashed line corre- 
sponds to configuration (b). 



VI. VISIBILITY OF THE INFORMATION 



The issue of the visibility of the information contained 
in the coherent SFG contribution, against the incoherent 
background, is a crucial one, especially for the observa- 
tion of the tails of the PDC X-shaped correlation. 



As already remarked, the expression (37) of the spec- 
tral coherence function of the SFG light, suggests that the 
visibility should be greatly enhanced by measuring the 
SFG light in the far-field of the second crystal, where the 
coherent contribution propagating in the forward direc- 
tion is separated from the incoherent background propa- 
gating over a broad angle. This issue will be discussed in 



Sec |VIB| We start here by considering the visibility of a 
measurement of the total number of SFG photons, which 
provides a useful estimation of the overall weight of the 
coherent component with respect to the incoherent one. 



case (a), Ay = in case (b)] is identical in the two 
configurations (see FigflQ^i). On the other hand, be- 
ing i/j meas (0, Ay, At) « in the Ay < region not 
too close to the central peak, relation ( [54] ) implies that 
Ilf G [0,Ay,At] « for Ay < and I™J[0, Ay, At] * 
4 x I^ h G [Ax = Ay, 0, A] for Ay > 0. The tails of the 
coherent SFG component, extending only in the Ay > 
region in configuration (b), are thus about four times 
more intense than in configuration (a). This behaviour is 
shown in FigjlO^c, which compares the temporal profile 
in the two configurations for the spatial shifts As = 50/im 
and As = — 50/mi respectively. 

We notice that configuration (b) offers the relevant 
benefit that the tails are more visible and can be there- 
fore detected more easily than in configuration (a). Al- 
though in principle configuration (b) does not allow the 
reconstruction of the symmetric X-shaped correlation as 
configuration (a), we found the following approximate 
empirical relation that links the results in the two con- 
figurations: 



I c s f G [Ax = Ay,0,At] 
1 



/§&[0, A V, At ] ~ V J sfg[0, - Ay, At] 



(55) 
for Ay ^ 0, 



this approximation holding as long as Ay is sufficiently 
large, i.e. far from the central peak. This relation should 
allow to infer from experimental data collected in con- 
figuration (b) the data that would be collected in con- 
figuration (a) , thus allowing the reconstruction of the 
X-shaped PDC correlation in configuration (b) . 



A. Bucket detection of the SFG photons 



The photon fluxes /gFG anc ^ ^sfg^A^j 0, At] given in 
Sec III , Eqs.(43), have been evaluated at the SFG crystal 



output face and are are uniform in the transverse plane 
because of the artifact of the plane-wave pump approxi- 
mation. Within this limit, the ratio 



V = 



r (coh) 
Z SFG 



[0,0, Az = 0] 



4 n G ) +4 C F O G ) [0,0,A, = 0] 



(56) 



defines therefore the visibility of the correlation peak at 
Ax = At = against the incoherent background either 
in the case of a near-field detection of the SFG intensity 
or assuming all the SFG photons are collected without 
discrimination (e.g. by using a bucket detector). 



Figure [TT] plots the visibility V as a function of the 
length of the SFG crystal V c for increasing PDC paramet- 
ric gain, starting from values of g below unity (FigfTTJi) 
up to very high values (FigfTTj}), the latter correspond- 
ing to the regime of the experiment reported in [4]. The 
contributions of b oth t he co herent and the incoherent 
components (43b) and (43a) have been estimated nu- 
merically with a Monte Carlo integration, assuming a 
550 nm FWHM temporal bandwidth is selected with a 
super-Gaussian frequency filter from 850 nm up to 1400 
nm. 



14 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



V 



0.01 1.0- 
0.1 

0.8 



^AAAAAAAAAAAAAZ 



(a) low-medium gain 



1000 2000 3000 4000 
4' (|um) 



0.3 0.6 
0.5 

0.4 
0.2 
0.0 



V 




(b) high gain 



1000 2000 3000 4000 

ic (m) 



^2000 



I™ h G [Ax = 0,At = 0] 



1000 



o 




1000 2000 3000 4000 



FIG. 11: Visibility of the coherent component with respect 
to the incoherent background as a function of the length of 
the SFG crystal l f c , for (a) medium-low and (b) high para- 
metric gain values (indicated in the figure). The length of the 
PDC crystal, l c = 4mm, is kept fixed. The selected temporal 
bandwidth extends from 850 nm up to 1400 mm. 



In a regime of low parametric gain, i.e. for g <C 1, the 
PDC spectral probability amplitude Fpdc(^) scales as 



FIG. 12: Estimation of the coherent and incoherent SFG pho- 



g [see Equation (14)], while the spectrum |y(u>)| 2 given 



For this reason and 



by Eq.(13) scales as g 2 

scale as g 2 and g 4 , respectively, and the visibility (|56|) is 



close to 100%. This behaviour can be physically under- 
stood by noticing that the probability of finding pair of 
photons that are not twins becomes exceedingly low in 
this regime. 

As the PDC parametric gain increases, the incoherent 
component rapidly increases and becomes comparable 
to the coherent one. As a consequence, V decreases to 
lower values, the degradation being particularly relevant 
in the l' c <C l c limit. For longer propagation distances 
in the SFG crystal, the generation of incoherent SFG 
photons becomes less efficient because of spatial walk-off 
and GVM, while the coherent component is not affected 
by those phenomena, as it will be further discussed in 
[T4] . The visibility remains therefore above 60% even for 
high PDC gains, as long as the propagation distance in- 
side the SFG crystal exceeds a few hundred micrometers. 
However, we notice that such visibility, although not neg- 
ligible, would make the detection of the full X-correlation 
very challenging in a practical implementation, because 
the tails would be hidden by the incoherent background. 



Fig 12 displays separately [0, 0, 0] and /gFG as a 



r (inc) 



of l' c , evaluated in the high-gain case (g 



function 

from Eq.(|43j). It shows that the coherent component in- 
creases almost quadratically up to l' c ~ 1 mm, while for 
larger propagation distance in the SFG crystal the co- 
herent up-conversion becomes less efficient. The initial 
quadratic behaviour reflects the result of Eq. (50) de- 
rived in the "thin SFG crystal" limit, where the filter- 
ing effect due to th e SF G spectral probability amplitude 
Fgp^(g, O) in Eq.(43b) is almost ineffective. 



By inspecting the expression of the typical bandwidths 

of Fgp^(g, and i^pDc(^i^) [ on ly the latter depends 
on the parametric gain according to Eq.(|20|)], we find that 



ton fluxes, /gFG an< ^ ^sfg J IP' 0' as a function of the SFG 
crystal length l' c . The PDC crystal parameters are g — 8 and 
l c 4000/im. 



r (coh) 



they become comparable for l' c = l c \Jn 2 /(ir 2 + g 2 ) « 
1.5 mm for g = 8, l c = 4 mm. For V c above this 
characteristic length the SFG spectral probability am- 
plitude .FsfcP b ecomes narrower than Fpdc(<L^), 
and the integral ( 43b ) spans a volume in the Fourier space 

oc /c~ 3 ^ 2 , so that does not scale any more as I' 2 . 

On the other side, the SFG incoherent component in- 
creases only linearly with l' c1 except for very short prop- 
agation distances for which approximation (51) holds. 



This is due to the fact that the up-conversion probability 
| ^sfg (^ w')\ 2 appearing inside the integral at r.h.s. of 
(43a) decays rapidly when w' ^ —w, because of the ef- 
fect of the spatial walk-off and the temporal group veloc- 
ity mismatch arising between the ordinary fundamental 
beam and the extraordinary up-converted beam. This 
issues will be further elucidated in a forthcoming related 
publication [T4] . 



B. Far-field detection of the SFG light 



According to the plane- wave pump result of (37), the 



spectral distribution of the PDC light consists of a coher- 
ent contribution, concentrated in the original plane-wave 
pump mode at q = and ft = 0, superimposed to an inco- 
herent background that spreads over a large bandwidth of 
spatio-temporal modes. The PWP approximations leads 
to artificial divergences (the ^-functions factors) that do 
not allow a direct evaluation of the visibility of the co- 
herent component against the incoherent background in 
the Fourier domain. For this reason, we implemented 
a full 3D+1 numerical simulation of our proposed setup 
(see Figj3|, and modeled the result of the detection of the 
SFG intensity distribution in the far field of the SFG crys- 
tal, as shown in fig|T3^. The generation of the broadband 
PDC field is simulated in the framework of the Wigner 
representation as described in [15], taking into account 
both the finite cross section w p and duration r p of the 
pump pulse, and the phase-matching conditions inside 
the two crystals (the full BBO Sellmeier dispersion rela- 



15 



tion [16 are used in the simulation). The propagation 
in the two crystals, described by Eqs.Q and ( |30| ), is 
simulated through a pseudo-spectral (split-step) method 
using a 256 x 256 x 512 numerical grid in the (x,t) — 
and the (g, t)— spaces. The parameters of the numeri- 
cal simulations are chosen to reproduce the conditions of 
the experiment being developed in Como [4] , which oper- 
ates in a pulsed regime of high parametric gain (g ~ 8). 
The broadband PDC field injected into the SFG crys- 
tal after the extraction of the pump beam undergoes the 
three- wave mixing process described by Eqs.(30) and the 



up-converted SFG field is mapped into the far-field with 



an f-f lens system, as shown schematically in Fig 13 1. We 
expect that a single stochastic realization of our simula- 
tions roughly reproduces the field distribution obtained 
from a single pump shot. 

The typical far-field intensity distribution of the up- 
converted SFG field obtained in the detection plane from 
a Gaussian pump pulse of waist w p — 600/im and du- 
ration r p = 1 ps is shown in Fig|l3j3,c. We verified that 
for the chosen PDC gain g=8 the injected PDC field (not 
shown in the figure) is only slightly depleted during prop- 
agation in the second crystal (< 0.1%), a feature which 
confirms the validity of the perturbative approximation 
d32bb for the SFG field used in the PWPA model. The 



narrow central peak results from coherent processes in 
which pairs of phase conjugate photons back-convert to 
a coherent field component reproducing partially the far- 
field distribution of the original pump beam. This peak 
is fixed, in the sense that it is reproduced identically in 
each stochastic simulation. The broad speckled back- 
ground instead is noisy and changes in each realization, 
giving rise on average to a broad distribution. It orig- 
inates from the incoherent processes where not phase- 
conjugated photon pairs are up-converted. These kind of 
simulations reproduce very closely the results on the far- 
field detection of SFG reported in [4] (see in particular 
Fig. 4 therein). 



The full scale plot in Fig[13|c shows that in such ideal 
conditions (perfect imaging, no dispersive optical ele- 
ments, no misalignments between the two crystals) the 
visibility of the coherent vs incoherent component for 
Ax = At = is very high, close to 100%. We can 
therefore expect that this far field detection scheme is 
well suited for performing the measurement of the full X 
correlation. 



Figure [14| simulates the reconstruction of the X-shaped 
correlation obtained by monitoring the number of pho- 
tons in the central peak of the SFG far-field as a function 
of the temporal delay At and the spatial shift Ax im- 
posed on the two twin PDC components (more precisely, 
it reports the number of SFG photons over a pixel 56 fim 
x 56 /im wide, in the focal plane of a lens with f=20 cm, 
placed at a focal distance from the SFG crystal). We 
verified that for the chosen pump pulse parameters, the 
retrieved structure in the (At, Ax)-plane is very close 
to that obtained by evaluating the coherent component 
( |47| ) within the PWPA (see e.g. Figj9^). Notice that 
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FIG. 13: (Color online) Numerical 3D+1 modeling of the 
experiment, (a) Detection scheme with a 2-f lens system to 
observe the SFG far field (f=20cm). (b) Far- field distribution 
of the SFG light for Ax = At = 0: the narrow central peak 
(truncated to 1% of its peak value) is the coherent component, 
arising from the up-conversion of phase conjugated photons. 
The incoherent contribution gives rise to the broad speckled 
distribution, (c) Full scale plot, showing that the visibility of 
the coherent component is close to 100%. g=8, w p = 600/im, 
r p = lps. The z-axis scale gives the estimated number of 
photons on the 56/mi x 56/xm pixels of the numerical grid. 



despite the fact that both the coherent and the incoher- 
ent components are displayed in Figfl4} the visibility of 
the tails of the correlation with respect to the speckled 
background is very high. The relative weight of the in- 
coherent background, about 0.2% of the coherent peak 
value, can be inferred from the height of the baseline of 
the temporal profile plotted in inset (d) for Ax = 0. 
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FIG. 14: (color online) Simulation of the reconstruction of 
the X-correlation via SFG under ideal conditions based on 
the far-field detection scheme of Fig 13 a. (a) Surface and (b) 



density plots of the number of SFG photons measured in the 
central pixel (at x = 0), obtained by scanning Ax and At on 
a 80 x 80 grid, (b) Temporal profiles for three different values 
of Ax showing the transition from the central peak to the 
double peak profile. The inset (d) shows the baseline of the 
Ax = profile, corresponding to the incoherent background 
with around 350 photons per pixel, about 0.2% of the peak 
value at Ax = At = 0. g = 8, w p = 600/xm, r p — lps. 



VII. FRAGILITY OF THE CORRELATION 
MEASUREMENT 

The previous simulations that display nearly 100% vis- 
bility for the correlation measurement have been ob- 
tained assuming a perfect imaging of the PDC source 
into the SFG crystal input face. The coherent component 
oc |^ meas (Ax, 0, At)| 2 is however strongly phase-sensitive 
and therefore extremely fragile against dispersion, imper- 
fect imaging conditions, as well as small misalignments 
between the two crystals. For this reason, the slightest 
imperfection in the imaging device that maps the PDC 
output plane into the SFG input plane deteriorates the 
coherent component of the SFG field which contains the 
correlation information. On the other hand, the incoher- 
ent component (43a) remains unaffected because of its 



phase-insensitive nature, so that the overall visibility of 
|?/; meas (Ax, 0, At) | 2 is strongly sensitive to those imper- 
fections. In the following we shall study the effects of the 
main sources of experimental imperfection. 



A. Temporal dispersion 

It is well known [7J [8j HQj |TT] that the scheme is very 
sensitive to the presence of any dispersive optical ele- 



ments. As an example, Fig. 15 shows the effect of a half 
millimeter thick slab of BK7 glass inserted in the prop- 
agation path between the PDC and the SFG crystals. 
It shows that not only the central correlation peak is 
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FIG. 15: Simulation of the effect of a dispersive element be- 
tween the PDC and SFG crystal obtained using the PWPA 
model. Spatio temporal correlation |V>meas| 2 - reconstructed 
(a) in ideal conditions, (b) in the presence of a 0.5 mm slab 
of BK7 glass between the PDC and the SFG crystals. 

broadened, but the X-shaped structure of correlation is 
strongly distorted. As a matter of fact, in order to com- 
pensate the detrimental effect of dispersion introduced by 
optical lenses, prisms have been used in the experimental 
works [3 [HI ED] ; a valid alternative approach, which we 
implemented in [4j[5] is rather to replace dispersive lenses 
by achromatic parabolic mirrors. 



B. Imperfect imaging 

Because of the smallness of the variation scale of 
the PDC X-correlation along the spatial dimension (< 
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10/im), we expect that the scheme shows very little toler- 
ance with respect to errors in the positioning of the SFG 
crystal with respect to the image plane of the telescopic 
system illustrated in Fig{l] Under ideal imaging condi- 
tions, the spatial width of the correlation peak is indeed 
determined by the collected PDC bandwidth considered 
in the simulation, i.e. 2g max w 0.6/im -1 , which gives a 
spatial correlation profile of width Axf w h m ~ 11 /im (see 
Fig|2}i and Figj^). 

Let us suppose for example that the second crystal 
entrance face is set at a distance / + Az[ mg from the 
second lens of the 4-f system. This error introduces an 
additional phase f actor to the transfer function product 
that enters in Eq.(43b), which is now given by 



ittAt iq x Ax 



-Az in 



(57) 

The condition for neglecting the diffraction term reads 



Az img <C £dof 



2 

max 



1 



n 2 



(58) 



where q mdiX is the maximal transverse wave- vector of the 
PDC emission in the collected bandwidth. For the tem- 
poral bandwidth 2^ max = 0.965 x 10 15 Hz assumed in the 
simulations, we have a tolerance on Az[ mg on the order 
of hundred micrometers. 

This is confirmed by our numerical modeling of the ex- 
periment. A first effect of an error in the imaging system 
is that the efficiency of the coherent up-conversion drops 
substantially with respect to the ideal configuration. Fig- 
ure [16] shows the far-field distribution of the SFG light 
(as it would be detected in the scheme of Fig 13 1), for in- 
creasing values of the error in the imaging plane Az[ mg . 
A displacement of 500 fim is sufficient to decrease the 
back-conversion efficiency of the twin photons by a fac- 
tor 12, and for Az mm = 2 mm the coherent peak becomes 
comparable in magnitude to the incoherent background 
(which is not sensitive to Azi mg ). 

Even more dramatical is the deterioration of the 
shape of the correlation function. Figure [TT] plots the 
spatio-temporal profile of the detected PDC correlation 
|^ meas (Ax, 0, At)\ 2 (row a and b) and of its temporal 
cross-section at Ax = (raw c) for increasing values 
of th e error Azi mg , evaluated from the PWPA result 
(43b). The density plots in row a) shows that a large 
broadening occurs along the spatial direction. This is 
not unexpected: with to respect an error in the imaging 
plane, the reconstructed spatial correlation as a func- 
tion of Ax behaves as the spatial resolution of an or- 
dinary optical image with respect to the depth of fo- 
cus of the imaging system, i.e. it roughly broadens as 

AXf w hm ^/l + A^img/^DOF" 

Less expected is perhaps the broadening of the correla- 
tion along the temporal direction, evidenced by raw c) in 
Fig 17 This temporal broadening is a clear consequence 
of the non-factorability of the twin-beam correlation in 
space and time. In turns, this is a consequence of the 
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FIG. 16: Simulation of the far- field distribution of the SFG 
light evaluated for increasing values of the error in the imaging 
plane positioning Azi mg . 



non-factorable character of phase matching: the phase 
matching condition [see Eq. (18)] 



A pdc (g, fi) « n 2 /n 2 D -q 2 /q 2 D = <-> k'+ti 2 = q 2 /h (59) 

can indeed be read as a compensation of the temporal 
dispersion experienced by twin photons inside the PDC 
crystal due to diffraction. Only if the entrance face of the 
second crystal is placed exactly in the image plane of the 
first crystal, this compensation occurs and the biphoton 
correlation is a nearly transform limited coherent sum of 
the phase-matched spectral modes. 

On the contrary, if the second crystal is misplaced with 
respect to the image plane, free propagation deterio- 
rates this exact compensation of dispersion and diffrac- 
tion. By making the simple assumption that only phase 
matc hed m odes contribute to the coherent SFG compo- 
nent (43b), we can make the substitution q 2 — >> Q 2 kik f { in 
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FIG. 17: (color online) (a) Behaviour of the reconstructed correlation | / mea s(Ax, 0, At)| 2 , normalized to its peak value, for 
increasing values of Azi mg , showing the rapid deterioration of the X-shape of correlation. The decrease of the peak height 
(normalized to the Azi mg = value) is shown in (b). The plot (c) shows the cross-section of the central peak along the 
temporal delay At, evidencing a broadening of the temporal correlation for increasing errors in the imaging plane 



the propagation term at the r.h.s. of (57). It transforms 

— I — 2 A^img 

then into the Vt— dependent phase factor e /w i 
which describes a quadratic dispersion-like chirp of the 
twin beams. 



C. Effect of misalignments of the two crystals 

The coherent component of the SFG light is also 
strongly sensitive to misalignments of the SFG crystal 
with respect to the PDC crystal orientation. We inves- 
tigated numerically the effect of a small tilt c)#o FG 7^ 
of the SFG crystal with respect to the PDC crystal ori- 
entation, the latter satisfying perfect phase-matching at 



degeneracy, i.e. 



0. Figure 
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shows how the 



peak of the coherent contribution , i.e [see Eqs.(45)-(46)]: 



rcoh 
i SFG 



[0,0,0] 



dw 



F PDC (w)F^\w) 



,(60) 



rapidly goes to zero as the tilt angle exceeds a few tenths 
of degrees. As expected, takes its maximum value 

when the two crystals are perfectly aligned, as Fpdc 



and i^sp^ , are peaked around the same phase-matching 
curves and their overlap integral is maximized. 

We can give a quantitative estimate of the angular tol- 
erance in the alignment of the two crystals by consider- 
ing explicitly the case when the second crystal is slightly 
tilted with respect to the first crystal, the latter being 
tuned for collinear phase matching. Denoting with # PDC 
(6 SFG ) the orientation angle of the PDC (SFG) crystals 
with the pump axis (i.e. the z-axis), and assuming the 
tilt angle between the two crystals 50 SFG = SFG - PDC 
is small, the following approximate relation between the 
pump mode wavenumbers in the two crystals holds 



uSFG 



ko - koPoSO 1 



SFG 



(61) 



where po = —dko / d qx \ UJ = 0JQ : q=o is the positive denned 
walk-off angle of the extraordinary wave at frequency 
ujq. Since by assumption A PDC = 2ki — ko =0, ac- 
cording to Eq.(61) a negative tilt S0 SFG < leads to a 
positive collinear phase-mismatch Aq FG = 2k\ — /cq FG = 
—kopoS0 SFG > in the second crystal. At degeneracy 
the phase-matched spatial modes in the second crystal 
lye therefore on a circumference of radius 



Qr 



WA SFG ^ 



WWo^ SFG |. (62) 
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FIG. 18: (color online) Bottom: peak of the coherent SFG 

intensity XgFC^IP' 0> ®\ as a function of the angular tilt A0q FG 
between the two crystals (normalized to the value it has for 
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-4.24, (b) S0q — 0° , A; 



V SFG;/ 



A§ FG Z' = 5.8. 



0, and (c) 



nSFG // 
"0 L c 



■ 0.07° 



component of the up-converted SFG field contains the 
desired information on the PDC correlation function. In 
particular, this coherent component can be expressed in 
the form of a convolution between the PDC biphotonic 
function in direct space, the quantity under investigation, 
with the corresponding probability amplitude describing 
the up-conversion process in the SFG crystal, the convo- 
lution being evaluated at the applied temporal delay and 
spatial shift between two conjugate PDC components. It 
is shown that the measured quantity retains the main 
features of the biphotonic correlation, namely its nonfac- 
torable X-shaped geometry and its strong localization in 
space and time. The latter, which can in principle re- 
duce to a few pump optical cycles pQ, is determined by 
the acceptance bandwidth of the up-conversion process. 
Finally, a fully 3D+1 numerical modeling of the optical 
setup that takes into account the pump pulse finite size 
has been implemented in order to provide a more real- 
istic simulation of the experiment being implemented in 
Como. This allowed us to show that, even in the high 
gain regime of PDC, the visibility of the correlation mea- 
surement con be close to 100% when the up-converted 
light is collected in the far-field of the SFG crystal, and 
for evaluating the tolerance of the phase-sensitive corre- 
lation measurement against imperfection of the imaging 
system. 



Close to degeneracy, the region of overlap of Fpdc an d 
can be expected to reduce drastically when qn 



(coh) 



SFG 



coincides with the first node of the PDC probability 
amplitude Fpdc(<T?^ = 0) along the spatial frequency 
which is found to be well approximate by qo 



V2(tt 2 +g 2 ) 1 / 4 q D [see Eq.(|20|)]. Following this criterium, 
the overlap integral (60) will be reduced by a large 



amount when the tilt angle exceeds the critical value 



56. 



SFG 



2v^ 



■r 



pokol' c 



(63) 



Taking e.g. g = 7 it gives us a tolerance 50^ FG = 0.2° in 
the alignment of the second BBO crystal with respect to 
the first. This value is in agreement with the numerical 
evaluation of [0, 0, 0] as a function of £#q FG shown 



in Fig 18 



VIII. CONCLUSIONS 

Our treatment shows that the SFG process represents 
a powerful tool for exploring the biphotonic correlation 
in the full spatio-temporal domain. Through a careful 
manipulation of both the spatial and the temporal de- 
grees of freedom of the PDC field emitted by the first 
^( 2 ) crystal, the proposed optical setup allows to retrieve 
the strongly localized X-shaped PDC correlation in the 
space time domain, predicted in [T]-[3] . The analytical re- 
sult obtained within the PWPA shows that the coherent 
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Appendix A: Derivation of Eqs.(26 )-(27) 



In this appendix we describe the effect of a rotating 
mirror set in the 2f-plane of a 4-f telescopic system. Let 
us consider a rotation by A<p of the mirror across the 
2/-axis (orthogonal to the figure plane) with respect to 
its normal position set at 45° with respect to the inci- 
dent beam (dashed line in Fig 19). Because of reflec- 



tion, the k- vector of an incident mode (<f, Q) with free 
space wavenumber kf(Q) — (uj\ + Q)/c and direction an- 
gle a = a,rcsm(q x / k f) in the (x, z)-plane is tilted by an 
additional angle 2A(j) with respect to the A(j) = direc- 
tion. Accordingly, the /c-vector transverse and longitu- 
dinal component in the plane orthogonal to the rotation 
axis, q x = kf sin a and k z = kf cos a, transforms accord- 
ing to the orthogonal transformation 



q' x = q x cos2A(/> + kf z sin2A(/> 



kfz cos2A0 — q x sin2A0 



(Ala) 
(Alb) 



20 




while that performed by the first lens can be written as 



if f f 
e H (g, CI) = -—61 \x = -— ft ) . 



(A5) 



Combining Eqs.( A3 )-( A5 ) together, we readily obtain the 
following transformation between the input and the out- 
put planes of the telescopic system 



FIG. 19: (color online) A rotation of A0 of the mirror with 
respect to its original position (dashed line) generates a tilt 
by 2A0 of the /c-vector with respect to the z-axis direction 
after reflection. 



ci(x, Q) = —bi(—x + Ax, ft) 



(A6) 



For small ang les, i.e. for q x /kf = sin a -C 1 and A<f> -C 1, 
relations ( Al ) reduces to 



q' x = q x + 2k f A<P, k' fz 



k 



fz 



(A2) 



where all powers of q x /k z and Aqb have been neglected 
except the linear terms. Under this approximation the 
relation between the reflected and the incident field, e\ r 
and en, can be written as 

eir(q x , q y , ft) = e u (q x - 2/c / (ft)A(/>, q y , ft) (A3) 

The field Fourier transformation performed by the second 
lens of the 4 — / system (see Figj3| can be written as 



c\ (x, ft) 



2tt/ 

ik f (n) 



ei, 



dx'e 1 "s x ' x ei r (x,Q) 



f 



-X, ^ 



(A4) 



where Ax = (2/A0, 0) denotes the transverse displace- 
ment of the field at the SFG crystal input plane produced 
by the mirror rotation. It is worth noticing that, under 
this approximation, the transverse shift is the same at all 
temporal frequencies. 



In the Fourier domain, Eq.(A6) reduces to the linear 
phase-shift transformation 



(A7) 



which coincides with input-output relation (25) affecting 
the q x > modes of beam +, except for the q — » —q 
reflection and the —1 factor introduced by the two lenses 
which has been omitted in the main treatment. 
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